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^ ■ Abstract. We present a detailed kinematical study and modeling of the emission 

■j-J ■ of the molecular cores at ambient velocities surrounding IRAS 21391+5802, an 

' intermediate-mass protostar embedded in IC 1396N. The high-density gas emis- 

sion is found in association with three dense cores associated with the YSOs 
• 1—1 . 

^ ■ BIMA 1, BIMA 2, and BIMA 3. The CS (J=5-^4) and CH3OH (J=5-i -^4_i) 

' 

' emission around BIMA 1 has been modeled by considering a spatially infinitely 

thin ring seen edge-on by the observer. From the model we find that CS is de- 
tected over a wider radii range than CH3OH. A bipolar outflow is detected in the 
CS (J=2^1) line centered near BIMA 1. This outflow could be powered by a yet 
undetected YSO, BIMA IW, or alternatively could be part of the BIMA 1 molec- 
ular outflow. The CS and CH3OH emission associated with the intermediate-mass 
protostar BIMA 2 is highly perturbed by the bipolar outflow even at cloud ve- 
locities, confirming that the protostar is in a very active stage of mass loss. For 
YSO BIMA 3 the lack of outflow and of clear evidence of infall suggests that both 
outflow and infall are weaker than in BIMA 1, and that BIMA 3 is probably a 
more evolved object. 



Key words. ISM: individual: IC 1396N, IRAS 21391-1-5802 — ISM: jets and 
outflows — ISM: molecules — radio lines: ISM — stars: formation 
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1. Introduction 

It is well known that Young Stellar Objects (YSOs) are embedded in dust and gas. For 
low-mass stars, the youngest and more embedded Class objects have most of their 
emission in large scale gas and dust structures, consistent with envelopes, while more 
evolved optically visible T Tauri stars have compact continuum emission, associated 
with disks, with little or no extended emission (e.g. Ohashi et al. 119911 119961 Chen et al. 
HM^l Looney et al. 000(1 . This well-defined evolutionary sequence found for low-mass 
stars, however, does not have a counterpart for higher mass stars, such as Herbig Ae 
and Be stars and their intermediate-mass precursors with masses in the range 2 Mq< 
Mi, < 10 Mq. For the intermediate- mass YSOs, the difficulty of identifying the youngest 
protostars, which are still in the infalling envelope phase, and the fact that they tend to 
be located at greater distances make them less studied and less well understood objects. 
The immediate vicinity of such protostars is a very complex environment, where the 
extended emission is usually resolved into more than one source when observed at high 
resolution (e.g. G173. 58+2.45: Shepherd & WatsonEOOS IRAS 21391+5802: Beltran et 
al. I2UU2L hereafter Paper I). In addition, the molecular outflows driven by intermediate- 
mass objects are more energetic. Thus, their interaction with the circumstellar gas and 
dust material surrounding the protostars is expected to be stronger and more dramatic, 
disrupting the envelopes and pushing away the dense gas at high velocities. 

For the intermediate-mass protostar IRAS 21391+5802 we conducted a detailed high 
angular resolution study by carrying out centimeter and millimeter continuum and spec- 
tral line BIMA observations (Paper I). The continuum emission at centimeter and mil- 
limeter wavelengths was resolved into three sources: BIMA 1, BIMA 2, and BIMA 3. 
The strongest source at millimeter wavelengths is BIMA 2, which is most likely the ob- 
ject associated with IRAS 21391+5802, and it has a circumstellar mass of ^ 5.1 M0. 
Regarding the CO emission, at least two molecular outflows were detected in the region, 
one of them is elongated along the north-south direction centered on BIMA 1, and the 
other one is a stronger east-west outflow centered at the position of the intermediate- 
mass protostar BIMA 2. The BIMA 2 CO outflow presents a complex structure and 
kinematics. While at high outflow velocities the outflow is clearly bipolar, at low out- 
flow velocities the blueshifted and redshifted emission are highly overlapping, and the 
strongest emission shows a V-shaped morphology. The CS and CII3OH emission from 
the BIMA 2 outflow exhibits two well differentiated and clumpy lobes, with two promi- 
nent northern blueshifted and redshifted clumps. The curved shape of the clumps and the 
spectral shape at these positions are consistent with shocked material. In addition, CS 
and CH3OH are strongly enhanced toward these positions with respect to typical quies- 
cent material abundances in other star-forming regions. These kinematical and chemical 

Send offprint requests to: Maria T. Beltran 
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evidences suggest that the clumps are tracing gas entrained within the surface of interac- 
tion between the molecular outflow and the dense ambient quiescent core, and that the 
morphology of the molecular outflow is a result of this interaction. As compared to the 
low-mass counterparts, the properties of the outflow driven by BIMA 2 are consistent 
with those of the outflows driven by low-mass young objects, and it fits well the corre- 
lations found for low-mass outflows by Cabrit & Bertout H1992;i . Anglada and 
Bontemps et al. In addition, the morphology and properties of the circumstel- 

lar structures around this intermediate-mass protostar are similar to those of low-mass 
counterpart Class objects. The other two sources in the region BIMA 1 and BIMA 3, 
have a mass of ^ 0.07 Mq each, and their small dust emissivity index is suggestive of 
grain growth in dense regions (e.g. Mannings & Emerson II 994|l . This fact together with 
the more compact appearance of their dust emission suggest that they could be more 
evolved low-mass objects. The results presented in Paper I analyzed and discussed the 
kinematics of the dense gas as it relates to the molecular outflows detected in the region. 
On the other hand, in order to continue our study of this intermediate-mass star-forming 
region, in this present work we analyze and model the kinematics of the gas at ambient 
velocities as it relates to the molecular cores in which BIMA 1, BIMA 2, and BIMA 3 
are embedded. The results of this study are presented here. 

2. Observations 

The details of the observations, carried out with the 10-antenna BIMA array^ between 
2000 March and 2001 January, are given in Paper I. The maps at 3.1 mm were done by 
using a robust weight of 0.5, with a resulting synthesized beam of 7'.'0 x 6'.'3, P.A.= —21°. 
Channel maps at 1.2 mm were also done with a robust weight of 0.5, with a resulting syn- 
thesized beam of 2'.'1 x I'.'S, P.A.= —3°, and the integrated intensity maps were done with 
a robust weight of 1, and the resulting synthesized beam is 2'.'7 x 2'.'4, P.A.= —16°. The 
spectral resolution of the observations was 0.30 km s^^ for CS ( J=2^1) at 97.981 GHz, 
and 0.12 km s^^ for CS (J=5^4) at 244.936 GHz and CH3OH (J==5_i ^4_i) at 
241.767 GHz. We note that for the velocity range at which the ambient gas arises (be- 
tween ~ —2 and 3 km s^^) the CS ( J=2^1) emission appears to be quite extended (see 
Fig. 1). This suggests that there is possibly a significant fraction of the total flux missed 
by the BIMA interferometer. The CS (J=2-^l) BIMA observations were done in the C 
array, so these observations should be sensitive to structures up to 60" (Wright 1996). The 
dense cores associated with BIMA 1, BIMA 2 and BIMA 3 have sizes (FWHM - 15") 
significantly below this value. The 1 mm observations (CH3OH and CS ( J=5^4)) were 
done using the B and C arrays, so they should be sensitive to structures up to 24". This 

^ The BIMA array is operated by the Berkeley-IUinois-Maryland Association with support 
from the National Science Foundation. 
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Fig. 1. Integrated intensity of the CS (J=2^1) emission over the velocity interval 
(—3,3) km s~^. The contom-s are 3, 5, 10, 15, 20, and 40 times 0.24 Jy beam^^ km s^^. 
The map has not been corrected for primary beam response. The crosses show the posi- 
tion of the 3.1 mm sources BIMA 1, BIMA 2, and BIMA 3 (Paper I), in order of right 
ascension. The dotted circumference represents the BIMA primary beam (50% attenua- 
tion level). The synthesized beam is shown in the lower left-hand corner. 

value is still larger than the CS ( J=2^1) sizes of the three dense cores. Therefore, the 3 
and 1 mm line emission from the three dense cores studied in this paper should not be 
significantly affected by missing flux. 

3. Results 

3.1. Lines observed 

CS (J=2^1), CS (J=5^4), and CH3OH (J=5_i ^4_i) were observed toward 
IRAS 21391-1-5802. Figure [T] shows the integrated emission of the CS ( J=2^1) at ambi- 
ent velocities (—3,3) km s^^. The strongest emission arises from a core associated with 
the YSO BIMA 2. West to this core, there is another core, which is associated with 
the YSO BIMA 1. Surrounding these two cores, there is extended and weaker emission 
spanning about 20-30" from the cores. The other YSO detected in the region, BIMA 3, 
is also associated with CS emission. Additionally, an elongation is also visible toward the 
northeast of BIMA 2. This emission is clearly seen in CS ( J=5— >4) as a clump located 
~ 15" from BIMA 2 (see Fig. 11 of Paper I). 
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Fig. 2. Spectra of the CS (J==2^1), CS (J=5->4) and CH3OH (J==5_i ^4_i) transi- 
tions toward the positions of BIMA 1, BIMA 2, and BIMA 3. The continuum emission 
has been subtracted. The 1-cr noise in 1 channel is 0.3 K for CS (J=2-^l), and 0.7 K for 
CS (J=5^4), and CH3OH (J==5_i ^4_i). The conversion factor is 2.92 K/Jybeam-i 
for CS (J=2^1), and 5.24 K/Jybeam-^ for CS (J=5-^4), and CH3OH (J=5_i ^4_i). 
The thick grey profiles are the Gaussian fits to the spectra of BIMA 1 and BIMA 3. 



CS ( J=5^4) and CH3OH ( J=5_i ^4_i) emission are also associated with the three 
YSOs BIMA 1, BIMA 2, and BIMA 3. The strongest and more extended emission for both 
species arises from the core associated with BIMA 2. Figure 2 shows the spectra obtained 
for CS (J=2^1), CS (J=5^4), and CH3OH (J=5_i ^4_i) toward the positions of 
BIMA 1, BIMA 2, and BIMA 3. 

3.2. Physical parameters 

Tabled lists the fitted line parameters for each transition for BIMA 1 and BIMA 3. The 
core associated with the YSO BIMA 2 is being strongly disturbed by the outflow even at 
velocities close to the systemic value, Vlsr — km s~^. As can be seen in the spectra in 
Figure 2, the BIMA 2 profiles are not Gaussian and show emission at high velocities. This 
is probably because the molecular outflow dominates the emission. In addition, both CS 
( J=2-^l) and ( J=5^4) lines also show a clear redshifted self- absorption feature. Thus, 
because of the difficulty to disentangle the cloud emission from the outflow emission, we 
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Table 1. CS (J=2^1), CS (J==5->4), and CH3OH (J==5_i -^4_i) lines and physical 
parameters obtained toward the positions of BIMA 1 and BIMA 3 around IRAS 21391- 
5802. 





Molecule & 








AV 


Ocs ^ ' 


71 r (h) 
A'/vir ^ ' 


Position 


Transition 


(K km s"-^) 


(K) 


( km s ) 


( km s ) 


(arcsec) 




BIMA 1 


CS (J=2^1) 


9.1 ±0.5 


6.9 ±0.5 


-0.54 ±0.03 


1.25 ±0.08 








CS (J=5^4) 


8.2 ±0.7 


5.8 ±0.7 


-0.74 ± 0.05 


1.33 ±0.13 


3.2 ±0.2 


0.31 ±0.12 




CH3OH 


4.8 ±0.8 


3.5 ±0.7 


-0.74 ±0.10 


1.27 ±0.27 






BIMA 3 


CS (J=2^1) 


8.1 ±0.5 


4.1 ±0.4 


±2.21 ±0.05 


1.86 ±0.15 








CS (J=5^4) 


9.7 ±0.7 


4.9 ±0.8 


±1.71 ±0.07 


1.87 ±0.17 


3.4 ±0.2 


0.3 - 4.6''=> 




CH3OH 


2.2 ±0.5 


2.1 ± 0.6 


±1.93 ±0.13 


0.98 ± 0.23s 







Deconvolved geometric mean of the major and minor axes of the 50% of the peak contour 
of the CS ( J=5~+4) map. 

Virial mass obtained from the CS (J=5-^4) using a AV of 0.5 km s~^, intrinsic line width 
(i.e. the line width corrected for kinematics) obtained from the collapsing ring model for BIMA 
1 (see §ETTJ. 

''^^ We cannot disentangle the kinematics of BIMA 3. Therefore we give the values for the virial 
mass assuming that the CS line width is the intrinsic line width (which gives an upper limit), 
and the value adopting AV = 0.5 km s~^, the intrinsic line width derived for BIMA 1. 

cannot study the properties of the molecular core associated with BIMA 2. The outflow 
clearly enhanced the integrated emission map shown in Figure ^ 

From the spectra (Fig. 2) and their fits (Table ^1 it is clear that while BIMA 1 has 
a systemic velocity of ^ —0.7 km s^^, the core associated with BIMA 3 is at a diS'erent 
velocity, Vlsr ~ 2 km s^^. In Tabled we also show the deconvolved size of the cores 
obtained from the CS (J—b^A) data, and the virial mass estimated from the intrinsic 
line width (i.e. the line width corrected for kinematics) obtained from the collapsing ring 
model for BIMA 1 (see ii l4.1.1|l . assuming a homogeneous spherical clump, and neglecting 
contributions from magnetic field and surface pressure. In such a case the virial mass can 
be computed from the expression (see e.g. MacLaren et al. I1988|l : 

V Mo y Vkpcy V arcsec / V km V 

where d is the distance, 750 pc for IRAS 21391+5802 (Matthews IWUjl . is the decon- 
volved size of the source, and AV is the intrinsic line width. It should be noted that the 
virial mass depends on the density profile. For a power-law density distribution of the 
type p oc r~P, the virial mass obtained from Eq. should be multiplied by a factor 
3(5 — 2p)/5(3 — p) which is < 1 for p < 3 (MacLaren et al. ll988|l . Thus, the values given 
in Tabled should be taken as upper limits. 
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Fig. 3. BIMA 1: Integrated intensity of the CS (J=2^1), CS (J=5^4), and CH3OH 
(J=5_i ^4_i) emission over the velocity interval (—1.8,0.1) km s^^ toward BIMA 1. 
The contours are 0.56 to 3.36 Jy beam^^ km s^^ in steps of 0.56 Jy beam^^ km s^^ 
for CS (J=2^1), -0.72, -0.48, -0.24, 0.24 to 1.44 Jy beam^i km s^^ in 
steps of 0.24 Jy beam^i km s^^ for CS (J=5^4), and -0.48, -0.24, 0.24 to 
0.96 Jy beam-i km s^^ in steps of 0.24 Jy beam^i km s^^ for CH3OH (J==5_i ^4_i). 
The cross shows the position of the 3.1 mm source BIMA 1 (Paper I). The synthesized 
beams are shown in the lower right-hand corner. The thick dashed-dotted line in the 
middle panel marks the axis of the north-south CO outflow driven by the YSO BIMA 1, 
at P. A. ~ 20° (Paper I). The thick solid hne marks the axis of the CS (J=2^1) blue 
and red lobes, with a P.A. ~ -5° (see S ET^ . 

4. Discussion 

4.1. BIMA 1 

As seen in the previous section the high-density gas emission at ambient velocities (—3, 3) 
km s~^ around IRAS 21391+5802 is found in association with three dense cores associated 
with the YSOs BIMA 1, BIMA 2, and BIMA 3. We conducted a kinematical study of 
each core that revealed important properties of the cores and of their evolutionary stage. 

Figure El shows the emission of CS (J=2->1), CS (J=5->4), and CH3OH 
(J=5_i ^4„i) toward BIMA 1, integrated for the velocity range (—1.8,0.1) km s~^. 
The core associated with the YSO BIMA 1, which is likely the powering source of the 
north-south bipolar molecular outflow detected in the region (Paper I), is well traced 
by the three transitions. On the one hand, CS ( J=5^4) and CH3OII peak at the same 
position as the 3.1 mm continuum emission (Paper I). On the other hand, as can be seen 
in this figure, the CS (J=2^1) emission peaks at a position shifted ^ 5" to the west 
from that of the CS (J=5^4) and CII3OH emissions. This shift in position is of the 
order of the CS (J=2— >1) beam (~ 7"), and thus, significant. Note that this is not an 
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Fig. 4. BIMA 1: {Bottom) PV plot of the CS (J=5^4) emission along the major axis, 
P.A.= 110° (left), and the minor axis, P.A.= 20° (right), of the BIMA 1 core. The position 
offset is relative to q(J2000) = 2l'^40'"39!89, (5(J2000) = 58°16'6'.'6. (Top) Synthetic 
emission of the PV plots along the major (left) and minor (right) axis of a ring with 
an infall velocity of 0.37 (i?/1000 AU)~"'^ km s^^, a rotation velocity of 0.15 km s^^, 
and inner and outer radius of 300 AU and 3750 AU. The contours are —0.3, 0.3 to 
1.05 Jy beam"^ (I'Sft) and to 1.20 Jy beam^^ (right) in steps of 0.15 Jy beam~^. 

instrumental or calibration effect, since the observations were carried out by configuring 
the BIMA digital correlator to observe simultaneously the 3.1 mm continuum emission 
and the CS (J=2->1) line. 

The morphology of the emission of the three species is also different. Despite the 
CS (J=2-^l) observations having a lower angular resolution than CS (J=5-^4), CS 
(J=:2^1) is tracing a more extended region, with a deconvolved size of ^ 14", about 
3 times larger than that traced by CS (J=5^4). All this suggests that CS (J==2^1) 
is tracing different material than CS (J=5— >4) and CII3OII, and thus, the physical 
properties and kinematics of the gas traced by them might be different. CS (J=5^4) 
and CH3OH would be likely tracing dense and hot gas associated with the YSO BIMA 1, 
while CS (J=2— >1) would be tracing more extended, less dense and cooler material. 

4.1.1. CS (J=5^4) and CH3OH (J=5_i ^4_i) 

As seen in Fig. 13 the integrated emission of CS (J=5^4) and CH3OII (J=5-i -^4_i) 
shows an elongated core at P.A.~ 110°, perpendicular to the outflow axis with P.A.~ 20°. 
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Fig. 5. BIMA 1: Same as Fig.gl for the CH3OH (J=5_i ^4_i) emission. The infalling 
ring has a radius of 1500 AU. The contom'S are —0.24, 0.24 to 0.60 Jy beam^^ (l-eft) and 
to 0.72 Jy beam^^ (right) in steps of 0.12 Jy beam~^. 



The deconvolved size of this core is (5'.'2±0'.'3) x (2'.'0±0'.'06), or (3900±230) x (1500±50) 
AU at the distance of the source. Schwartz et al. H1991II have mapped in C^^O ( J=1^0) 
a similar structure toward the position of BIMA 1. If this morphology corresponds to a 
disk-like structure seen nearly edge-on, the ratio between the major and minor axes of 
the emission gives an inclination angle of ^ 70°. 

The kinematics of the gas surrounding BIMA 1 can be seen in the position-velocity 
(PV) cuts done along the major axis (110°) and the minor axis (20°) of the dense core 
(see Figs. 0] and IS)). For both molecules, the PV plot along the major axis shows a 
structure roughly symmetric both in velocity, respect to a velocity of —0.7 km s^^, and 
in position, respect to the center. The plot along the minor axis shows two blobs at the 
central position, symmetric respect to a velocity of —0.7 km s^^ as well. 

For CH3OH (J=5-i ^4_i), the PV plot along the major axis shows a ringlike mor- 
phology, suggesting that the gas is collapsing (or expanding) radially and that the emis- 
sion arises from a ringlike or toroidal structure. The CS (J=5— +4) PV plot, however, 
shows a different morphology. The highest velocities (relative to —0.7 km s^^) of the 
core appear only near the central position, while for the largest position offsets, ~ ±4", 
the relative velocities are smaller. This can be interpreted as the emission of superpos- 
ing rings of different radius, with an infall (or expansion) velocity that decreases with 
increasing radius. 
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We considered a model of a spatially infinitely thin disk with infall and rotation, 
similar to that of Ohashi et al. ifPU^ for L1554 and that of Girart et al. PM|l for HH 
SON. Regarding the inclination, we considered a disk seen edge-on by the observer. We 
think that the inclination i of the disk is small because the PV signature seen in the 
data agrees with i ~ 0° (see Figs. 0] and Isj: for both molecules the emission along the 
major axis extends roughly from —5" to +5", while along the minor axis it is compact, 
and for CH3OH the PV plot along the major axis shows a characteristic ring shape. In 
case the disk is not edge-on, one should apply a correction factor proportional to sin^ i. 
We assumed that the infall velocity was a power law of the radius, with index —0.5 (the 
value expected for free fall), V^^f oc i?^"-^. Since the fit showed that the rotation velocity 
was very small, we assumed it to be independent of radius. In addition, we assumed 
that the intensity was proportional to R^^ . In case of optically thin emission, this would 
correspond to a disk with a surface density proportional to R^^ . The free parameters of 
the model were the inner and outer radii of the disk, i?inn and i?out, the infall velocity at 
a fixed radius (taken to be 1000 AU), V^nf, and the rotation velocity, V^ot- The synthetic 
PV diagrams were computed along the projected major and minor axes of the disk, and 
convolved with a Gaussian in the PV plane, with a resolution of 2" in position and 0.5 
km s"^ in velocity, to compare with the observed diagrams for both molecules. 

The observed and synthetic PV diagrams are shown in Figs. 0] and |S1 . We explored 
the parameter space of the model and found that it was possible to found a good fit 
to the four observed PV diagrams simultaneously, for the same values of the velocities 
(^nf = 0.37 km s"^, T^ot < 0.15 km s~^), and different values of the inner and outer 
radii of the disk for each molecule (CS: i?inn = 300 AU, i?out = 3750 AU; CH3OH: 

i?inn = i?out = 1500 AU). 

Assuming that the infalling material is in free fall, the dynamical mass required is 
given by Afdyn — ^nf-^/^- From the value and uncertainty of the infall velocity obtained 
from the fit, we can estimate it to be Mdyn = 0.16 ± 0.05 M0. The virial mass estimated 
from the CS ( J=5^4) intrinsic line width is Afvir — 0.31 Mq (Table^) for a homogeneous 
density distribution. However, for typical density distributions, p oc r^P, with p — 2.0- 
1.5, the virial mass is 0.18-0.24 Mq. Thus, taking into account the uncertainties of the 
model, such values are consistent with the dynamical mass. Regarding the circumstellar 
mass, the value estimated from the dust emission (Paper I) is 0.07 Mq, roughly half 
of the dynamical mass. If the observed velocities are infall, this would require a total 
(circumstellar plus protostar) mass similar to the dynamical mass, and thus, the protostar 
should have a mass similar to the circumstellar mass. This is consistent with BIMA 1 
being a late Class or Class I protostar (Andre et al. l2()()()|l . Thus, the data available are 
consistent with the BIMA 1 core being in collapse. 

As mentioned above the radii of the ring were i?inn = 300 AU and -Rout = 3750 AU 
for CS, and i?inn = -Rout = 1500 AU for CH3OH. The value of the outer radius can be 
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affected by sensitivity. Its value is constrained by the extension of the PV map in the 
position direction. However, this should not affect significantly the estimate of the mass 
of the disk because the power-law index for the density obtained is —1.0. Note that the 
CS emission spreads over a larger range of radii, up to 3750 AU, while that of CH3OH 
is detected only at a radius of 1500 AU. This probably traces molecular abundance 
variations with radius (CH3OII being confined within a narrow range of radii from the 
protostar), or different excitation conditions for both species. In the latter case the non 
detection of CS and CH3OH at radii smaller than i?inn could be attributed to the fact 
that higher energy transitions than those observed are excited for the higher temperature 
near the protostar. Note that a CH3OII knot visible at a velocity of ~ —1.5 km s~^ 
(FiglSJ is not fitted by the model. The velocity of this knot does not fit the kinematic 
structure of the circumstellar gas. However, the knot may be part of the blueshifted gas, 
since the blueshifted lobe of the outflow encompasses the position of BIMA 1 (see next 
section and Fig. EJ. 

4.1.2. CS (J=2^1) 

CS {J—2—>-l) spectrum in Fig. 2 shows marginally wings of emission, especially at 
blueshifted velocities with respect to the Vlsr of the core, —0.7 km s^^. Figure 
shows the CS (J=2— >1) emission integrated for the velocity intervals corresponding to 
the dense core, to the blueshifted and to the redshifted emission. The blueshifted and 
redshifted lobes are located roughly symmetrically on both sides of the peak of the CS 
(J=2^1) core emission. This suggests that the CS ( J=2^1) emission traces outflowing 
gas at moderate velocities. The axis of this outflow has a P.A.~ —5°. The PV diagram 
along this axis is shown in the right panel of Fig. 1^1 where the blueshifted wing is clearly 
visible toward the north of the core, and the redshifted wing toward the south. The spec- 
tra of the CS ( J=2^1) emission taken at the center of each lobe (see top panel of Fig.EJ 
shows the presence of blueshifted and redshifted material, confirming the existence of an 
outflow traced by the CS (,/=2^1) emission. Assuming that the CS emission is in LTE 
and is optically thin, we calculated a mass of roughly 0.02 Mq for the gas associated 
with the outflow, for a excitation temperature of 10 K, and a CS abundance of SxlO^'^ 
(e.g., NGC 1333 IRAS2A: Jorgensen et al. 12004(1 . The momentum and kinetic energy of 
the BIMA IW outflow is 0.05 Mq km s~^ and 1 x 10^^ ergs, respectively, which are 
clearly lower than those of the BIMA 1 and BIMA 2 outflows traced by CO (Paper I). 

The center of symmetry of the CS (J=2^1) outflow is located ^ 9" west of BIMA 1, a 
shift that is slightly larger than the synthesized beam diameter along the right ascension. 
An inspection of the channel maps at the ambient gas velocities shows that the peak 
emission is located between BIMA 1 and the center of symmetry of the CS (J=2^1) 
outflow. This is also shown in the map of the integrated emission for the ambient gas 
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Fig. 6. BIMA 1: {Bottom left) Overlay of the integrated intensity of the CS 
(J=2— >1) emission for the velocity interval (—1.5,0.8) km s~^ (core emission, 
greyscale), for (—3,-1.5) km s~^ (blueshifted emission, thick solid contours), and for 
(0.8,3) km s~^ (redshifted emission, thick dashed contours). Greyscale levels range 
from to 1.9 Jy beam~^ km s~^. Contours levels are 2, 3, 4, 5, 6, and 7 times 
0.15 Jy beam~^ km s~^. The white cross shows the 3.1 mm position of BIMA 1 (Paper 
I). The dashed-dotted line marks the axis of the north-south CO outflow driven by the 
YSO BIMA 1, at P. A. ~ 20° (Paper I). The thin sohd hne marks the axis of the CS 
(J=2^1) blue and red lobes, with a P.A. ~ -5°. {Bottom right) CS (J=2-»l) PV dia- 
gram along the solid line of left panel. The origin of the PV diagram is at the position 
a(J2000) = 2l'^40"39!00, ^(J2000) = 58°16'5'.'l. Contour levels are -2, 2, 3, 4, 5, 6, 7, 
and 9 times 0.2 Jy beam~^. {Top) Spectra of the CS (J=2— >1) emission at the centers 
of the blue and red lobes. 

velocities (Fig. 6). All this suggests that the CS (J=2— >1) outflow could be powered by 
a yet undetected YSO, BIMA IW, located at the center of symmetry of this outflow 
(see Fig. 6 caption for the absolute coordinates). Based on the millimeter continuum 
observations (Paper I), the 3a upper limit for the dust mass of this embedded object 
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Fig. 7. BIMA 2: Overlay of the integrated intensity of the CS (J=2^1) emission (left), 
and the CS ( J=5^4) emission (right) over the velocity interval (—3, 3) km s^^ [greyscale) , 
and the velocity field of the dense gas [contours) around BIMA 2. The solid contours 
depict blueshifted velocities with respect to the systemic velocity, and the dashed contours 
depict redshifted velocities. The contours are —1.5, —1.2, —0.9, —0.6, and —0.3 km s^^ for 
blueshifted velocities [solid contours), and 0.3, 0.6, 0.9, 1.2, and 1.5 km s^^ for redshifted 
velocities [dashed contours). Greyscale levels are 0.2 to 10.2 Jy beam^^ km s^^ in steps 
of 2 Jy beam~^ km s~^. The cross shows the 3.1 mm position of the YSO BIMA 2. The 
solid line marks the axis of the east-west CO outflow driven by the YSO, at P.A. ~ 60° 
(Paper I). 



would be ~ 0.04 Mq. Note that if BIMA IW has associated a compact core traced by 
CS (J=2->1), as the other YSOs in this region do, BIMA 1 and BIMA IW could not 
be spatially resolved because of the CS (J=2^1) angular resolution, and the blended 
emission will appear to peak between the two sources, as it is the case. 

Alternatively, we cannot discard that the CS ( J=2^1) outflow is not an independent 
outflow but it is part of the outflow associated with BIMA 1 and detected in CO ( J=l— >0) 
at higher velocities (see Paper I). The northern and southern lobes of this CO outflow 
are redshifted and blueshifted, respectively, which is the reverse situation than that of 
the CS outflow. Yet, this could be explained if the BIMA 1 outflow axis is near to the 
plane of the sky. In any case, in this scenario the CS (J=2— >1) bipolar emission would 
possibly be tracing the denser part of the low velocity cavity of the BIMA 1 outflow. 

It is clear that further high angular resolution (^ 3") and high sensitivity observations 
at millimeter and submillimeter wavelengths would be needed to confirm the presence of 
the BIMA IW embedded young source and its association with the CS ( J=2— >1) outfiow. 
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Fig. 8. BIMA 2: PV plot of CS (J=2^1), CS and CH3OH (J=5_i ^4_i) 

along the axis of the molecular outflow driven by the YSO BIMA 2, P.A.~ 60°, with 
the offset position relative to the 3.1 mm position of the YSO BIMA 2 (Paper I). The 
contours are —0.8, 0.8 to 2.8 Jy beam^^ in steps of 0.4 Jy beam^^ for CS (J=2^1), 
-0.3, 0.3 to 1.8 Jy beam^i in steps of 0.15 Jy beam^i for CS (J=5^4), and -0.3, 0.3 
to 0.9 Jy beam-i in steps of 0.15 Jy beam^ifor CH3OH (J=5_i ^4_i). 



4.2. BIMA 2 

BIMA 2 is the intermediate-mass protostar associated with IRAS 21391+5802 and is 
powering a strong east- west molecular outflow (Paper I). As seen in Fig. ^ the molecular 
core associated with this YSO is the strongest core in the region. The CS and CH3OH 
emission are being strongly disturbed by the outflow even at ambient velocities. Thus, it 
is not possible to study in much detail the kinematics of the core because the presence 
of outflowing gas along the line-of-sight makes any interpretation of the motions very 
difficult. This can be seen in Fig. [3 where the first-order moment (intensity weighted 
mean Vlsr) of the CS ( J=2^1) and ( J=5^4) emission around BIMA 2 shows a velocity 
gradient across the source in the direction of the molecular outflow (P.A.~ 60°). This 
velocity is correlated with the molecular outflow; that is, blueshifted velocities toward the 
southwest and redshifted velocities toward the northeast. This can also be seen in the PV 
diagrams of Fig. Hthat show CS (J=2^1), CS (J=5->4), and CH3OH (J=5_i ^4_i) 
along the axis of the molecular outflow. 

For this core, the CS (J=2^1) and CS (J=5— >4) integrated emission peaks close to 
the position of 3.1 mm continuum source (Paper I). However, Fig.|Slshows that around 
BIMA 2 the CS ( J=2^1) emission peaks at - -0.5 km s"^ whereas the CS ( J=5^4) 
emission peaks at ~0.5 km s^^. This along with the clear self-absorption feature seen in 
both transitions at the position of the dust source (Fig. 2) is likely due to a high optical 
depth of ambient gas emission. Therefore, this suggests that most of the CS emission 
at ambient velocities is not tracing the dense core around BIMA 2 but the interaction 
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Fig. 9. BIMA 2: Integrated intensity of the CH3OH (J=5_i ^4_i) emission toward 
BIMA 2. The cross shows the 3.1 mm position of the YSO BIMA 2. Greyscale levels are 
to 6 Jy beam^^ km s^^ in steps of 1 Jy beam^^ km s^^. Solid line marks the axis 
of the east- west CO outflow driven by the YSO BIMA 2, at P. A. ~ 60° (Paper I). The 
synthesized beam is shown in the lower left-hand corner. 



layer of the molecular outflow with the dense core. This is supported by the kinematical 
correlation of the CS emission at ambient velocities with the molecular outflow. 

The CH3OH (J=5_i ^4_i) emission peak is offset by ~ 2" with respect to the CS 
(J=5-^4) emission peak, of the order of the beam size, in the outflow direction (see 
Fig. El). The PV cuts of Fig. |S1 show a deflcit of CII3OII emission toward the position of 
the 3.1 mm source for all velocities. This molecule is usually produced by low- velocity 
shocks induced by the molecular outflow and it is more often found in association with 
the outflowing gas (e.g. Bachiller et al. I1998|l . Thus, this suggests that the chemistry of 
CH3OH is probably being affected by shocks between the molecular outflow driven by 
BIMA 2 and the dense core at ambient velocities. However, the spatial offset between 
the CS and CII3OII suggests that they trace different outflowing material, due either to 
different chemical conditions or to different excitation conditions. 

In summary, for this intermediate-mass protostar, both CS and CII3OH arc not good 
tracers of the dense core at ambient velocities but of the interaction of the molecu- 
lar outflow with the core. This scenario is different from the low-mass case, where the 
high-density gas surrounding the protostar is clearly tracing the dense core at ambi- 
ent velocities. Thus, in order to study the kinematics of the dense core associated with 
BIMA 2, observations of different high-density tracers are needed. 




a(2000) 0((2000) 

Fig. 10. BIMA 3: Integrated intensity of the CS ( J=2^1) and CS ( J=5^4) emission 
over the velocity interval (0.5, 3.0) km s^^ toward BIMA 3. The contours are 0.56 to 
2.24 Jy beam^^ km s^^ in steps of 0.56 Jy beam^^ km s^^ for CS (J=2^1), and 
-0.48, -0.24, 0.24 to 1.92 Jy bcam-i km s^^ in steps of 0.24 Jy beam^i km s^^ for CS 
(J=5^4). The cross shows the position of the 3.1 mm source BIMA 3 (Paper I). The 
synthesized beams are shown in the lower left-hand corner. 

4.3. BIMA 3 

The molecular core associated with the YSO BIMA 3 has been clearly detected in CS 
(J=2^1) and CS (J=5^4), and marginally, in CH3OH (J=5_i ^4_i) (see Figs. 2 and 
I10() . However, no molecular outflow was found in association with it in Paper I. This 
core is well defined but it is clearly at a velocity different to that of the BIMA 1 and 
BIMA 2 cores, as shown by the spectra (Fig. 2) and by the different moment maps of the 
velocity field. The Vlsr of this core is around ^ 2 km s^^. The CS (J=5^4) integrated 
emission shows a northwest-southeast fiattened core, which has a deconvolved size of 
4'.'1 ± O'.'l X 2'.'9 ± 0'.'04, or 3100 ± 75 x 2200 ± 30 AU at the distance of the source, at a 
P.A.2± 150° (see Fig-UHl). 

CS (J=5^4) PV cuts done roughly in the direction of this elongation (P.A.~ 150°) 
and in the direction perpendicular (P.A.~ 60°) are shown in Figure [TTl These diagrams 
are similar to those of BIMA 1 (Fig.0}, and thus, they could be indicative of collapse or 
expansion of the molecular core. The CS ( J=5^4) PV cuts could not be properly fitted 
by our simple model (the CH3OH emission is too marginal to include it in the fit). This 
suggests a non-uniform density distribution of the circumstellar material around BIMA 
3, which might be due to fragmentation and the presence of subcores within BIMA 3. 
Alternatively, the infall process in BIMA 3 could be weaker than in BIMA 1. 
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Fig. 11. BIMA 3: PV plot of the CS (J=5->4) emission along the major axis, P.A.= 150° 
(left) , and the minor axis, P.A.= 60°, (right) of the BIMA 3 core. The position offset is 
relative to a(J2000) = 2li^40'"42!91, (5(J2000) = 58°16'00'.'3. The contours are -0.3, 0.3 
to 1.2 Jy beam^^ in steps of 0.15 Jy beam^^. 

4.4. Evolutionary stage of the dense cores 

A detailed study of the emission of the molecular cores at ambient velocities surrounding 
the intermediate-mass protostar IRAS 21391+5802 has allowed us to infer or confirm the 
evolutionary stage of the cores studied in Paper I. 

For YSO BIMA 1, we suggest a model for the CS and CII3OH emission, proposing 
that it is a low-mass object undergoing collapse. From the CO ( J=1^0) data (Paper I) 
we estimated that the dynamical timescalc for the north-south molecular outflow centered 
at the position of BIMA 1 is of the order of 10^ years. Such a value is consistent with 
dynamical timescales estimated for outflows driven by Class objects (e.g. Andre et 
al. HUM Richer et al. IMIHjl . Near BIMA 1 the detection of a CS (J=2-^l) bipolar 
outflow suggests the existence of a yet undetected YSO, BIMA IW. The dynamical 
timescale estimated for the CS (J=2^1) outflow is 1-2x10^ years, which suggests that 
the powering object is also a very YSO. However, further observations are required in 
order to confirm this new source and derive its properties. 

We found that the CS and CII3OII emission associated with YSO BIMA 2 is strongly 
disturbed by the outflow, indicating that the protostar is in a very active stage of mass 
loss, which confirms that this intermediate-mass object is a very young stellar object, 
probably a Class object. The dynamical timescale estimated for the CO ( J=1^0) east- 
west molecular outflow driven by BIMA 2 (Paper I) is of the order of 10^ years, which 
also confirms that BIMA 2 is a very YSO. 

Regarding YSO BIMA 3, it is much less massive (M ~ 0.07 M0; Paper I) than 
BIMA 2 (A/ ~ 5.1 Mq; Paper I), and its emission is more compact. This core is associated 
with maser emission (Tofani et al. I1995|l and centimeter continuum emission (Paper I). 
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Both are clear signposts of star formation and indicate that the outflowing phase has 
already started. Furthermore, the detection of CH3OH, although marginal, indicates 
the occurrence of shocks in this core, likely produced by an outflow, although can also 
indicate mantle evaporation due to the heating from the protostar. However, although 
those observations suggest that BIMA 3 is a YSO, no large scale molecular outflow 
emission has been detected, and no clear and conclusive infall evidence is visible from 
the PV plots (see Fig. 111(1 . Thus, both infall and outflow are probably weaker for BIMA 3 
than for BIMA 1, which suggests that BIMA 3 would be a more evolved object, possibly 
a low-mass Class I object. The molecular outflows powered by such objects are usually 
much weaker than the ones driven by Class objects (e.g. Bontemps et al. 11996(1 . which 
would explain the non detection of the molecular outflow toward BIMA 3. 

None of the YSOs BIMA 1, BIMA 2, or BIMA 3 has been detected at J, H, or 
Kg band through the infrared 2MASS survey. However, in the immediate vicinity of the 
YSOs there are other infrared (proto)stellar sources detected at Kg band (2.17 fim). This 
means that in a relatively small region there are very young stellar objects, detected at 
millimeter wavelengths (the BIMA sources) together with more evolved young objects 
(the 2MASS sources). 

5. Conclusions 

We conducted a kinematical study and modeled the gas emission of the molecular cores 
at ambient velocities surrounding IRAS 21391+5802, an intermediate-mass protostar 
embedded in the core of IC 1396N. The high-density gas has been found in association 
with three dense cores associated with the YSOs BIMA 1, BIMA 2, and BIMA 3. 

For the YSO BIMA 1, the integrated emission of CS (J=5->4) and CH3OH 
(J=5_i ^4_i) shows an elongated core at a P.A. of ~ 110°, perpendicular to the out- 
flow axis. We modeled the gas emission by considering a spatially infinitely thin ring with 
infall and rotation seen edge-on by the observer, with the intensity and infall velocity 
distributions being power laws of the radius. The best fit was obtained with the inten- 
sity proportional to -R^^, an infall velocity Vinf — 0.37 (-R/1000 AU)^"'^ km s^^, and a 
rotation velocity l^ot < 0.15 km s^^. Assuming that the infalling material is in free fall, 
the dynamical mass required is 0.16 Mq, which is similar to the total (circumstellar + 
protostar) mass in the core. 

The CS ( J=2^1) emission toward BIMA 1 peaks at a position ~ 5" west from that 
of the 3.1 mm continuum, CS (J=5^4), and CH3OH emissions. The high-velocity CS 
(J=:2^1) emission is tracing outflowing gas at moderate velocities (~ 2 km s^^). The 
center of symmetry of this CS outflow is near the CS (J=2^1) emission peak. This 
outflow could be powered by a yet undetected YSO, BIMA IW, or alternatively could 
be part of the BIMA 1 molecular outflow. 
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The CS and CH3OH omission associated with the intcrmcdiatc-mass protostar 
BIMA 2 arc not tracing the dense core but they are highly perturbed by the bipolar 
outflow even at velocities close to the systemic value. This confirms that the protostar 
is in a very active stage of mass loss, and also supports the idea that intermediate-mass 
outflows interact more dramatically with the dense gas surroimding the protostar than 
those driven by low-mass protostars. It also suggests that CS and CH3OH are likely not 
good tracers for studying the properties of dense cores associated with intermediate-mass 
protostars. 

The core associated with BIMA 3 is at a systemic velocity different to that of the 
BIMA 1 and BIMA 2 cores. The CS (J=5^4) integrated emission shows a northwest- 
southeast flattened core. The lack of molecular outflow and of clear evidence of infall 
suggests that both outflow and infall are weaker than in BIMA 1. Thus, BIMA 3 is 
probably a more evolved. Class I, object. 

CS (J=5-^4) PV cuts are similar to those of BIMA 1, and thus, they could also be 
indicative of collapse of the molecular core, with an infall velocity of ~ 0.4 km at 
a radius of 1000 AU. This infall velocity corresponds to a dynamical mass of 0.14 Mq, 
comparable to the total mass. 
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